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ABSTRACT The regulation of cellular function is often controlled at the level 
of gene transcription. Such genetic regulation usually consists of interacting 
networks, whereby gene products from a single network can act to control their 
own expression or the production of protein in another network. Engineered 
control of cellular function through the design and manipulation of such networks 
lies within the constraints of current technology. Here we develop a model 
describing the regulation of gene expression, and elucidate the effects of noise 
on the formulation. We consider a single network derived from bacteriophage 
A, and construct a two-parameter deterministic model describing the temporal 
evolution of the concentration of A repressor protein. Bistability in the steady- 
state protein concentration arises naturally, and we show how the bistable regime 
is enhanced with the addition of the first operator site in the promotor region. 
We then show how additive and multiplicative external noise can be used to 
regulate expression. In the additive case, we demonstrate the utility of such 
control through the construction of a protein switch, whereby protein production 
is turned "on" and "off" using short noise pulses. In the multiplicative case, we 
show that small deviations in the transcription rate can lead to large fluctuations 
in the production of protein, and describe how these fluctuations can be used 
to amplify protein production significantly. These novel results suggest that 
an external noise source could be used as a switch and/or amplifier for gene 
expression. Such a development could have important implications for gene 
therapy. 
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Introduction 



Regulated gene expression is the process through which cells control fundamental functions, 
such as the production of enzymatic and structural proteins, and the time sequence of this 
production during development ||1], |2|. Many of these regulatory processes take place at the 
level of gene transcription |j| , and there is evidence that the underlying reactions governing 
transcription can be affected by external influences from the environment [|]]. 

As experimental techniques are increasingly capable of providing reliable data pertaining 
to gene regulation, theoretical models are becoming important in the understanding and 
manipulation of such processes. The most common theoretical approach is to model the 
interactions of elements in a regulatory network as biochemical reactions. Given such a 
set of chemical reactions, the individual jump processes (i.e., the creation or destruction 
of a given reaction species) and their associated probabilities are considered. In its most 
general form, this often leads to a type of Monte Carlo simulation of the interaction proba- 
bilities ||. Although this approach suffers from a lack of analytic tractability, its strength 
is its completeness - fluctuations in species' concentrations are embedded in the modeling 
process. These internal fluctuations are important for systems containing modest numbers 
of elements, or when the volume is small. 

Rate equations originate as a first approximation to such a general approach, whereby 
internal fluctuations are ignored. These deterministic differential equations describe the evo- 
lution of the mean value of some property of the set of reactions, typically the concentrations 
of the various elements involved. The existence of positive or negative feedback in a regu- 
latory network is thought to be common ||, and, within the reaction framework, feedback 
leads to nonlinear rate equations J7|. 

Noise in the form of random fluctuations arises in these systems in one of two ways. As 
discussed above, internal noise is inherent in the biochemical reactions. Its magnitude is 
proportional to the inverse of the system size, and its origin is often thermal. On the other 
hand, external noise originates in the random variation of one or more of the externally set 
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control parameters, such as the rate constants associated with a given set of reactions. If the 
noise source is small, its effect can often be incorporated post hoc into the rate equations. In 
the case of internal noise, this is done in an attempt to recapture the lost information embod- 
ied in the rate equation approximation. But in the case of external noise, one often wishes 
to introduce some new phenomenon where the details of the effect are not precisely known. 
In either case, the governing rate equations are augmented with additive or multiplicative 
stochastic terms. These terms, viewed as a random perturbation to the deterministic pic- 
ture, can induce various effects, most notably the switching between potential attractors 
(i.e., fixed points, limit cycles, chaotic attractors) ||. 

While impressive progress has been made in genome sequencing and the understand- 
ing of certain qualitative features of gene expression, there have been comparatively few 
advancements in the quantitative understanding of genetic networks. This is due to the 
inherent complexity of such biological systems. In this work, we adopt an engineering ap- 
proach in studying a solitary gene network. We envision that a plasmid, or genetic applet ||, 
containing a small, self-contained gene regulatory network, can be designed and studied in 
isolation. Such an approach has two distinct advantages. First, since the approach is inher- 
ently reductionist, it can make gene network problems tractable and thus more amenable 
to a mathematical formulation. Secondly, such an approach could form the basis for new 
techniques in the regulation of in vivo gene networks, whereby a genetic applet is designed 
to control cellular function. 

In this paper, we develop a model describing the dynamics of protein concentration in 
such a genetic applet, and demonstrate how external noise can be used to control the net- 
work. Although our results are general for networks designed with positive autoregulation, 
we ground the discussion by considering an applet derived from the promotor region of bac- 
teriophage A. Since the range of potentially interesting behavior is wide, we focus primarily 
on the steady-state mean value of the concentration of the A repressor protein. This choice 
is motivated by experiment; detailed dynamical information is still rather difficult to obtain, 
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as are statistical data concerning higher moments. We show how an additive noise term 
can be introduced to our model, and how the subsequent Langevin equation is analyzed by 
way of transforming to an equation describing the evolution of a probability function. We 
then obtain the steady-state mean repressor concentration by solving this equation in the 
long time limit, and discuss its relationship to the magnitude of the external perturbation. 
This leads to a potentially useful application, whereby one utilizes the noise to construct a 
genetic switch. We then consider noise at the level of transcription, where noise enters the 
formulation in a multiplicative manner. As in the additive case, we transform to an equation 
describing a probability distribution, and solve for the steady-state mean concentration as a 
function of noise strength. Finally, we demonstrate how such a noise source can be used to 
amplify the repressor concentration by several orders of magnitude. 



A Model for Repressor Expression 

In the context of the lysis- lysogeny pathway in the A virus, the autoregulation of A repressor 
expression is well-characterized |10| |. In this section, we present two models describing the 



regulation of such a network. We envision that our system is a plasmid consisting of the 
Pji — Prm operator region and components necessary for transcription, translation, and 
degradation. 

Although the full promotor region in A phage contains the three operator sites known 
as OR1, OR2, and OR3, we first consider a mutant system whereby the operator site OR1 
is absent from the region. The basic dynamical properties of this network, along with 
a categorization of the biochemical reactions, are as follows flO|| . The gene cl expresses 



repressor (CI), which in turn dimerizes and binds to the DNA as a transcription factor. In 
the mutant system, this binding can take place at one of the two binding sites OR2 or OR3. 
(Here, we ignore nonspecific binding.) Binding at OR2 enhances transcription, which takes 
place downstream of OR3, while binding at OR3 represses transcription, effectively turning 
off production. 



The chemical reactions describing the network are naturally divided into two categories 
- fast and slow. The fast reactions have rate constants of order seconds, and are therefore 
assumed to be in equilibrium with respect to the slow reactions, which are described by rates 
of order minutes. If we let X, X2, and D denote the repressor, repressor dimer, and DNA 
promoter site, respectively, then we may write the equilibrium reactions 

m 
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D + X 2 




DX 2 


D + X 2 




DX* 








DX 2 + X 2 




DX 2 X 2 



where the DX 2 and DX 2 complexes denote binding to the OR2 or OR3 sites, respectively, 
DX 2 X 2 denotes binding to both sites, and the Ki are forward equilibrium constants. We let 
K 3 = G\K 2 and K4, = cr 2 K 2 , so that o\ and a 2 represent binding strengths relative to the 
dimer-OR2 strength. 

The slow reactions are transcription and degradation, 

DX 2 + P DX 2 + P + nX (2) 
X H A 

where P denotes the concentration of RNA polymerase and n is the number of proteins per 
mRNA transcript. These reactions are considered irreversible. 

If we consider an in vitro system with high copy-number plasmidsQ, we may define con- 
centrations as our dynamical variables. Letting x = [X], y = [X 2 ], d = [D], u = [DX 2 ], 
v = [DX%], and z = [DX 2 X 2 ], we can write a rate equation describing the evolution of the 
concentration of repressor, 

x = —2kix 2 + 2k-iy + n ktPo u ~ kdX + t (3) 



*This assumption is necessary since the number of relevant molecules per cell is small in vivo. Since there 
are many cells, we could alternatively use state probabilities as dynamical variables describing an in vivo 
system. 



where we assume that the concentration of RNA polymerase po remains constant during 
time. The parameter r is the basal rate of production of CI, i.e., the expression rate of the 
cl gene in the absence of a transcription factor. 

We next eliminate y, u, and d from Eq. ([3]) as follows. We utilize the fact that the 
reactions in Eq. (jl]) are fast compared to expression and degradation, and write algebraic 
expressions 

y = K lX 2 (4) 

u = K 2 dy = K 1 K 2 dx 2 

v = 0\K 2 dy = u\K\K 2 dx 2 

z = a 2 K 2 uy = a 2 (K 1 K 2 ) 2 dx 4 

Further, the total concentration of DNA promoter sites dx is constant, so that 

d T = d + u + v + z = d(l + (1 + a^K^x 2 + a 2 K 2 K 2 x 4 ) (5) 

Under these assumptions, Eq. @ becomes 

. nk t p Q d t K l K 2 x 2 

X ~ 1 + (1 + o x )K x K 2 xi + a 2 K 2 K 2 x* + T W 

Without loss of generality, we may eliminate two of the parameters in Eq. (0) by rescaling 
the repressor concentration x and time. To this end, we define the dimensionless variables 



x = x\/K\K 2 and t = t{r\/K\K 2 ). Upon substitution into Eq. (]3|), we obtain 



2 

CX X /rr\ 

x = -, r - s + 1 (7) 

1 + 1 + (Jik 2 + a 2 x 4 



where the time derivative is with respect to t and we have suppressed the overbar on 
x. The dimensionless parameter a = nk t p dT/r is effectively a measure of the degree in 
which the transcription rate is increased above the basal rate by repressor binding, and 



7 = kd/ ir\/K\K 2 ) is proportional to the relative strengths of the degradation and basal 
rates. 
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For the mutant operator region of A phage, we have o\ ~ 1 and o 2 ~ 5 [[T^, |Tl|], so that 
the two parameters a and 7 in Eq. (0) determine the steady-state concentration of repressor. 
For this equation, there are two types of behavior. For one set of parameter values, we have 
monostability, whereby all initial concentrations evolve to the same fixed-point value. For 
another set, we have three fixed points, and the initial concentration will determine which 
steady state is selected. Additionally, in the multiple fixed-point regime, stability analysis 
indicates that the middle fixed point x m is unstable, so that all initial values x < x m will 
evolve to the lower fixed point, while those satisfying x > x m will evolve to the upper. This 
bistability arises as a consequence of the competition between the production of x along with 
dimerization and its degradation. For certain parameter values, the initial concentration is 
irrelevant, but for those that more closely balance production and loss, the final concentration 
is determined by the initial value. 

Graphically, we can see how bistability arises in Eq. (0) by setting a x 2 / (1 + 2 x 2 + 5 x 4 ) = 
7 x — 1. In Fig. 1A we plot the functions a x 2 /(l + 2 x 2 + 5 x 4 ) and 7 x — 1 for fixed a 
and several values of the slope 7. We see that for 7 small (whereby degradation is minimal 
compared with production), there is one possible steady-state value of x (and therefore CI). 
As we increase 7 above some critical value 7^, we observe that three fixed-point values 
appear. As we increase 7 still further beyond a second critical value 7^, the concentration 
"jumps" to a lower value and the system returns to a state of monostability. 

The preceding ideas lead to a plausible method whereby the system may be experimen- 
tally probed for bistability. We envision that a is fixed by the transcription rate and DNA 
binding site concentration, and that the degradation parameter 7 is an adjustable control. 
Beginning with a low initial value of 7 = 70 = 5, we slowly increase the degradation rate. 
The effect is illustrated in Fig. IB. We see that as 7 is slowly increased, the concentration 
of CI slowly decreases as the system tracks the fixed point. Then, at the moment when 7 
is greater than 7^, the concentration abruptly jumps to a lower value, followed by a further 
slow increase. Now suppose we reverse course, and begin to decrease 7. Then the system 
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will track along the lower fixed point until a point when 7 is greater than 7^. At this point, 
the system will again jump, this time to a higher fixed-point value. The trademark of hys- 
terisis is that the two jumps, one when increasing 7 and the other when decreasing, occur 
for different values of 7. 

As is well-known, the full operator region of A phage contains three sites. We turn briefly 
to the effect of the additional site OR1 on the above network. In order to incorporate its 
effect, Eq. ([[]) must be generalized to account for additional equilibrium reactions. This 
generalization amounts to the incorporation of dimer binding to OR1 [[L(J , and permutations 
of multiple binding possibilities at the three operator sites. Then, using known relationships 
between the cooperative binding rates, the above steps can be repeated and an equation 
analogous to Eq. (|7|) constructed. We obtain 

a (2 x 2 + 50 x 4 ) 

X ~ 25 + 29 x 2 + 52 x 4 + 4 x 6 ~ 7 X ( ) 

As can be seen, the addition of OR1 has the effect of changing the first term on the right- 
hand side of the equation. While this augmentation does not affect the qualitative features 
of the above discussion, one important quantitative difference is depicted in Fig. IB. In this 
figure, we see that the addition of OR1 has a large effect on the bistability region, increasing 
the overall size of the region by roughly an order of magnitude. Additionally, the model 
predicts that, while the drop in the concentration of repressor at the first bifurcation point 
will be approximately the same in both cases, the jump to the higher concentration will be 
around five times greater in the system containing OR1. Finally, since one effect of a larger 
bistable region is to make the switching mechanism more robust to noise, these results are 
of notable significance in the context of the lysogeny-to-lysis switching of A phage. 

Additive Noise 

We now focus on parameter values leading to bistability, and consider how an additive ex- 
ternal noise source affects the production of repressor. Physically, we take the dynamical 
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variable x described above to represent the repressor concentration within a colony of cells, 
and consider the noise to act on many copies of this colony. In the absence of noise, each 
colony will evolve identically to one of the two fixed points, as discussed above. The pres- 
ence of a noise source will at times modify this simple behavior, whereby colony-to-colony 
fluctuations can induce novel behavior. 

An additive noise source alters the "background" repressor production. As an example, 
consider the effect of a randomly varying external field on the biochemical reactions. The 
field could, in principle, impact the individual reaction rates |12|, and since the rate 
equations are probabilistic in origin, its influence enters statistically. We posit that such an 
effect will be small and can be treated as a random perturbation to our existing treatment; 
we envision that events induced will affect the basal production rate, and that this will 
translate to a rapidly varying background repressor production. In order to introduce this 
effect, we generalize the aforementioned model such that random fluctuations enter Eq. (§) 
linearly, 

x = f(x) + £(t) (9) 

where f(x) is the right-hand side of Eq. (H), and is a rapidly fluctuating random term 
with zero mean (< >= 0). In order to encapsulate the rapid random fluctuations, we 
make the standard requirement that the autocorrelation be "5-correlated" , i.e., the statistics 
of are such that < >= D5(t — t'), with D proportional to the strength of the 

perturbation. 

Eq. (||) can be rewritten as 

* = +{W (10) 

where we introduce the potential 4>(x), which is simply the integral of the right-hand side of 
Eq. (0). <j)(x) can be viewed as an "energy landscape", whereby x is considered the position 
of a particle moving in the landscape. One such landscape is plotted in Fig. 2A. Note that 
the stable fixed values of repressor concentration correspond to the minima of the potential 
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<f) in Fig. 2A, and the effect of the additive noise term is to cause random kicks to the particle 
(system state point) lying in one of these minima. On occasion, a sequence of kicks may 
enable the particle to escape a local minimum and reside in a new valley. 

To solve Eq. ([Top, we introduce the probability distribution P(x,t), which is effectively 
the probability of finding the system in a state with concentration x at time t. Given Eq. . 



a Fokker-Planck (FP) equation for P(x,t) can be constructed ||T4| 

d t P(x,t) = -d x (f(x)P(x,t)) + ^d x 2 P(x,t) (11) 

We focus here on the value of the steady-state mean (ssm) concentration. To this end, we 
first solve for the steady-state distribution, obtaining 

P s (x) = Ae-^ (x) (12) 

where A is a normalization constant determined by requiring the integral of P s (x) over all x 
be unity. In Fig. 2B, we plot P s (x), corresponding to the landscape of Fig. 2 A, for two values 
of the noise strength D. It can be seen that for the smaller noise value the probability is 
distributed around the lower concentration of repressor, while for the larger noise value the 
probability is split and distributed around both concentrations. This is consistent with our 
conceptual picture of the landscape: low noise will enable only transitions from the upper 
state to the lower state as random kicks are not sufficient to climb the steep barrier from the 
lower state, while high noise induces transitions between both of the states. Additionally, 
the larger noise value leads to a spreading of the distribution, as expected. 

Using the steady-state distribution, the steady-state mean (ssm) value of 
given by 

POO 

<x> ss = / xAe~^ m dx (13) 



Jo 

In Fig. 2C, we plot the ssm concentration as a function of D, obtained by numerically 
integrating Eq. ( |T3"D and transforming from the dimensionless variable x to repressor concen- 
tration. It can be seen that the ssm concentration increases with D, corresponding to the 
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increasing likelihood of populating the upper state, as discussed previously with respect to 
Figs. 2A and B. 

Figure 2C indicates that the external noise can be used to control the ssm concentration. 
As a candidate application, consider the following protein switch. Given parameter values 
leading to the landscape of Fig. 2A, we begin the switch in the "off" position by tuning the 
noise strength to a very low value. This will cause a high population in the lower state, and 
a correspondingly low value of the concentration. Then at some time later, consider pulsing 
the system by increasing the noise to some large value for a short period of time, followed 
by a decrease back to the original low value. The pulse will cause the upper state to become 
populated, corresponding to a concentration increase and a flipping of the switch to the "on" 
position. As the pulse quickly subsides, the upper state remains populated as the noise is 
not of sufficient strength to drive the system across either barrier (on relevant time scales). 
To return the switch to the off position, the upper-state population needs to be decreased to 
a low value. This can be achieved by applying a second noise pulse of intermediate strength. 
This intermediate value is chosen large enough so as to enhance transitions to the lower 
state, but small enough as to remain prohibitive to upper-state transitions. 

Figure 2D depicts the time evolution of the switching process for noise pulses of strengths 
D = 1.0 and D = 0.05. Initially, the concentration begins at a level of [CI] = 10 nM, 
corresponding to a low noise value of D — 0.01. After six hours in this low state, a 30- 
minute noise pulse of strength D = 1.0 is used to drive the concentration to a value of 
[CI] ~ 58 nM. Following this burst, the noise is returned to its original value. At 11 hours, 
a second 90-minute pulse of strength D = 0.05 is used to return the concentration to its 
original value. 

Multiplicative Noise 

We now consider the effect of a noise source which alters the transcription rate. Although 
transcription is represented by a single biochemical reaction, it is actually a complex sequence 

11 



of reactions |T(J, and it is natural to assume that this part of the gene regulatory sequence 
is likely to be affected by fluctuations of many internal or external parameters. We vary the 
transcription rate by allowing the parameter a in Eq. (H) to vary stochastically, i.e., we set 
a — > a + £(£). In this manner, we obtain an equation describing the evolution of the protein 
concentration x 

x = h{x) + £{t)g{x) (14) 

where h(x) is the right-hand side of Eq. @, and 

M _ 2 x 2 + 50 x 4 

9{X> 25 + 29 x 2 + 52 x 4 + 4 x 6 1 ' 

Thus, in this case, the noise is multiplicative, as opposed to additive, as in the previous case. 

Qualitatively, we can use the bifurcation plot of Fig. 3A to anticipate one effect of allowing 
the parameter a to fluctuate. Such a bifurcation plot is yet another way of depicting the 
behavior seen in Fig. 1A; it can be seen that for certain values of a there is one unique 
steady-state value of repressor concentration, and that for other values there are three. 
To incorporate fluctuations, if we envision a to stochastically vary in the bistable region 
of Fig. 3A, we notice that the steep top branch implies the corresponding fluctuations in 
repressor concentration will be quite large. This is contrasted with the flat lower branch, 
where modest fluctuations in a will induce small variations. In order to verify this observation 
quantitatively, we simulated Eq. flT4]), the results of which are presented in Fig. 3B. Beginning 
with repressor concentration equal to its upper value of approximately 500 nM, we notice 
that the immediate fluctuations are quite large even though a varies by only a few percent 
(Fig. 3A). Then, at around 700 minutes, the concentration quickly drops to its lower value, 
indicating that the fluctuations envisioned in Fig. 3A were sufficient to drive the repressor 
concentration to the dotted line of Fig. 3A and off the upper branch (across the unstable 
fixed point). The final state is then one of very small variation, as anticipated. 

As in the previous section, the steady-state probability distribution is obtained by trans- 
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forming Eq. fllil) to an equivalent Fokker-Planck equation JT4 



d t P(x, t) = -d x (h(x) + ^g(x)g'(x))P(x, t) + ^d 2 u g 2 (u)P(x, t) (16) 

where the prime denotes the derivative of g(x) with respect to x. We again solve for the 
steady-state distribution, obtaining 

P s {x) = Be~^ m{x) (17) 

As before, the steady-state distribution can be used to obtain the ssm concentration. 

Although not originating from a deterministic equation like that of Eq. ([7]), the function 
4> m (u) in Eq. (|D]) can still be viewed as a potential. We now consider parameter values 
leading to one such landscape in Fig. 3C. This landscape implies that we will have two 
steady-state repressor concentrations of approximately 5 and 1200 nM. This large difference 
is due to the largeness of the parameter a, implying that repressor "induced" transcription 
amplifies the basal rate by a large amount. (Since dr enters in the numerator of the definition 
of a, one could construct such a system experimentally with a high copy-number plasmid). 
This feature suggests that multiplicative noise could be used to amplify protein production, 
as described in the following example. We begin with zero protein concentration and very 
low noise strength D, leading to a highly populated lower state and low overall concentration. 
Then, at some later time, we pulse the system by increasing D for some short interval. This 
will cause the upper state to become quickly populated as it is easy to escape the shallow 
valley of the landscape and move into the large basin. In Fig. 3D, we plot the temporal 



evolution of the mean repressor concentration obtained from the simulation of Eq. fllil) . We 
see that the short noise pulse at around 20 hours indeed causes the concentration to increase 
abruptly by over three orders of magnitude, making this type of amplification an interesting 
case for experimental exploration. 
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Discussion 



From an engineering perspective, the control of cellular function through the design and 
manipulation of genetic regulatory networks is an intriguing possibility. In this paper, we 
have shown how external noise can be used to control the dynamics of a regulatory network, 
and how such control can be practically utilized in the design of a genetic switch and/or 
amplifier. Although the main focus of this work was on a network derived from the promotor 
region of A phage, our approach is generally applicable to any autoregulatory network where 
a protein-multimer acts transcription factor. 

An important element of our control scheme is bistability. This implies that a neces- 
sary criterion in the design of a noise-controlled applet be that the network is poised in a 
bistable region. This could potentially be achieved by methods such as the utilization of a 
temperature-dependent repressor protein, DNA titration, SSRA tagging, or pH control. 

Physically, the noise might be generated using an external field. Importantly, it has 
been claimed that electromagnetic fields can exert biological effects ||15|| . In addition, recent 
theoretical [13 and experimental |12|] work suggests a possible mechanism whereby an electric 
field can alter an enzyme-catalyzed reaction. These findings suggest that, although there is 
global charge neutrality, an external field can interact with local dipoles which arise through 
transient conformational changes or in membrane transport. 

Current gene therapy techniques are limited in that transfected genes are typically ei- 
ther in an "on" or "off" state. However, for the effective treatment of many diseases, the 
expression of a transfected gene needs to be regulated in some systematic fashion. Thus, the 
development of externally-controllable noise-based switches and amplifiers for gene expres- 
sion could have significant clinical implications. 
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Figure Captions 

FIG. 1. Bifurcation plots for the variable x and concentration of A repressor. (A) Graph- 
ical depiction of the fixed points of Eq. (^), generated by setting x = and plotting 
a x 2 /(l + 2x 2 + 5a; 4 ) and the line jx — 1. As the slope 7 is increased, the system traverses 
through a region of multistability and returns to a state of monostability. (B) Hysterisis 
loops for the mutant and nonmutant systems obtained by setting x = in Eqs. ([?[) and (|8|). 
Beginning with concentrations of 35 nM for the mutant system and 85 nM for the nonmutant 
system, we steadily increase the degradation parameter 7. In both systems, the concentration 
of repressor slowly decreases until a bifurcation point. In the mutant (nonmutant) system, 
the repressor concentration abruptly drops to a lower value at 7 ~ 16 (7 ~ 24). Then, 
upon reversing course and decreasing 7, the repressor concentration increases slowly until 7 
encounters a second bifurcation point at 7 ~ 14 (7 ~ 6), whereby the concentration immedi- 
ately jumps to a value of 15 nM (mutant) or 70 nM (nonmutant). The subsequent hysterisis 
loop is approximately 10 times larger in the nonmutant case. Parameter values are a = 50, 
Ki = 0.05 nM" 1 , and K 2 = 0.026 nM" 1 for the mutant system, and K 2 = 0.033 nM" 1 for 



the nonmutant system [10 



FIG. 2. Results for additive noise with parameter values a = 10 and 7 = 5.5. (A) The energy 
landscape. Stable equilibrium concentration values of Eq. (^) correspond to the valleys at 
[CI] = 10 and 200 nM, with an unstable value at [CI] = 99 nM. (B) Steady-state probability 
distributions for noise strengths of D = 0.04 (solid line) and D = 0.4 (dotted line). (C) The 
steady-state equilibrium protein concentration plotted versus noise strength. The concentra- 
tion increases as the noise causes the upper state of (A) to become increasingly populated. 
(D) Simulation of Eq. (Q) demonstrating the utilization of external noise for protein switch- 
ing. Initially, the concentration begins at a level of [CI] = 10 nM corresponding to a low 
noise value of D = 0.01. After six hours, a large 30-minute noise pulse of strength D = 1.0 
is used to drive the concentration to 58 nM. Following this pulse, the noise is returned to 
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its original value. At 11 hours, a smaller 90-minute noise pulse of strength D = 0.04 is used 
to return the concentration to near its original value. The simulation technique is that of 
Ref. HI7|. 



FIG. 3. Results for multiplicative noise. (A) Bifurcation plot for the repressor concentration 
versus the model parameter a. The steep upper branch implies that modest fluctuations in 
a will cause large fluctuations around the upper fixed value of repressor concentration, while 
the flat lower branch implies small fluctuations about the lower value. (B) The evolution of 
the repressor concentration in a single colony, obtained by simulation of Eq. fll^l) . Relatively 
small random variations of the parameter a (~ 6%) induce large fluctuations in the steady- 
state concentration until around 700 minutes and small fluctuations thereafter. (C) Energy 
landscape for parameter values a = 100 and 7 = 8.5. (D) Large-scale amplification of the 
protein concentration obtained by simulation of Eq. (|14]) . At 20 hours, a 60-minute noise 
pulse of strength D — 1.0 is used to quickly increase the protein concentration by over three 
orders of magnitude. The parameter values are the same as those in (C). 
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Figure 1 - Hasty et al. 




Figure 2 - Hasty et al. 




Figure 3 - Hasty et al. 



